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Abstract 

In 1987, J. H. Elton m, has proved the first fundamental result in convergence of IFS, the El¬ 
ton’s Ergodic Theorem. In this work we prove the natural extension of this theorem to the projected 
Hutchinson measure fj, a associated to a GIFSpdp <S = (X, (< f>j : X m —> X)j = o,i,...,n-i, {Pj)j= o,i,...,n-i), 
in a compact metric space (X, d). More precisely, the average along of the trajectories x„(a) of the 
GIFS, starting in any initial points xo, £ X satisfies, for any / € C(X, R), 

i JV_1 r 

lim a? X 

Jv ^+°° N 1 Jx 

for almost all a £ fl = {0,1,..., n— 1} N , the symbolic space. Additionally, we give some examples and 
applications to Chaos Games and Nonautonomous Dynamical Systems defined by finite difference 
equations. 


Keywords: Generalized iterated function system with probabilities, Markov operator, Hutchinson mea¬ 
sure, Ergodic Theorem, Iterated Function Systems, Dynamical Systems, Chaos Games. 
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Introduction 

In 2008, A. Mihail and R. Miculescu [16], has introduced the Generalized Iterated Function Systems 
(GIFS, for short). They prove that there exist a fractal attractor and give estimates of the rate of 
convergence for contractive GIFS. In 2009, Alexandru Mihail HE. has considered the Hutchinson measure 
associated to a GIFS with place dependent probabilities (GIFSpdp for short) that generalizes the classical 
Hutchinson measure, the invariant measure, associated to an Iterated Function System (IFS for short). 

In this work, the central idea is to extend the GIFS to an IFS on a bigger space, that we call Extended 
GIFS. From this IFS we get an Extended Hutchinson measure. In the rest of the paper, we set up the 
properties of the extended Hutchinson measure and his relation with the original GIFS. 

Using the extension, we prove an ergodic theorem that extends the classic Elton’s ergodic theorem 
for IFS, Elton EH, to the Hutchinson measure associated to a GIFS. From our results, we prove a Chaos 
game theorem for GIFS. As an application, we get some results on the stability and the asymptotic 
behavior of nonautonomous dynamical systems defined by finite difference equations. Also, we show haw 
to recover some properties of Gibbs measures for Holder potentials through an appropriate GIFS that is 
builded from an expansive endomorphism. 

The paper has three sections. In the Section |T| we recall the basic facts about GIFS and GIFSpdp, 
that we will use in the rest of the paper. In the Section [2] we introduce the extension of a GIFS and prove 
the Ergodic Theorem (Theorem [24]). Section [3] is devoted to applications. The main goal here is to prove 
the Chaos Game Theorem (Theorem |2U1) , allowing us to draw the attractor of the extended GIFS and 
his projection. 

We believe that the tools we develop will be very useful to forthcoming works and for other researchers 
in this area. The ergodic theorem we prove, and his consequences, represents a real advance in the 
understand of GIFS. 


1 Background on GIFS and GIFSpdp 

In this section, we will recall the basic definitions and results on the theory of GIFS. See R. Miculescu m , 
for more details and the notation. We notice that the word “generalized” has been used in several different 
ways representing more general components of a classical IFS. Here, generalized, means that functions 
are from X m to X instead X to X. 


1.1 Generalized iteration function system (GIFS) 

Let (X, d) be a compact El metric space (typically X = [0,1], X = {0,1,..., d— 1} N , X = [0,1] N , etc). 
Consider the topology on X m given by 

doo ((xi,x 2 , Xm ), ( 2 / 1 ,2/2,-,2/rn)) = . max d{xi,yi), 

1 = 1,...,771 

then ( X m , doo) is also a compact metric space. 

Definition 1. A (continuous) generalized iterated function system (GIFS) of degree m is a (finite) 
family S of continuous functions tpj : X m —>■ X, denoted S = (X, o...n-i)- 

See Secelean [T8HT9] for the analogous theory for countable GIFS. In order to avoid technicalities 
we assume that m — n — 2 that is, two maps in X 1 2 (see Remark [6] for additional details). 
So the standard GIFS is S — (X, (</>j)j=o,i) where </>o? <fii '• X 2 —» X. We recall that, 

Lip(X 2 , X) = {/ | d(f(x 0 ,y 0 ), f(x 1 ,y 1 )) < Cd{{x 0 ,y 0 ), (xi,yi)), C := Lip{f)} 


and 

Lipa^X 2 , X) = {f | d(f(x 0 ,yo), f(xi,yi)) < ad(x 0 ,x i) + bd(y 0 , yi), a, b > 0}. 
From now on we will assume the contraction hypothesis: 


1 We assume compactness to avoid technicalities. Many of the results that we present here are true if (X,d) is just 

complete. It is sufficient because the measures we use are always supported in the attractors that are compact sets. 
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El - Each (f>j : X 2 — > X is in Lip aj ^ j (X 2 , X) and aj + bj < 1. In particular, all the tpj are Lipschitz 
contractions and Lip{cj)j) = aj + bj. 

As usual, we denote /C(X) C V(X) = 2 X , the family of compact subsets of X. Moreover, Af*(X) = 
/C(X)\{0}. 

Definition 2. Given f : X 2 — > X we define the associated set function Ff : V*(X) 2 —> V*(X) by 
Ff(K\,Kfi) = f(Ki x X2). Also, we define the function F$ : X*(X) 2 — > X*(X) associated to S by 
Fs(K\, K2) = U;=o 1 ^4 >d (Xi, Kf)- A set Y C X is self-similar (or fractal) with respect to S if Fg(Y, Y ) = 

Y. 

The map Fg is sometimes called, Fractal operator, Barnsley’s Function or Hutchinson’s operator, in 

the literature. From Mihail and Miculescu [16], Theorem 3.5, we know that under the hypothesis El, there 

exists a unique attractor A(<S) £ /C*(X) for the GIFS that depends continuously on <pj. That is, A(5) 

is self-similar (Fg(A(S), A(S)) = A(S)) and, for any H 0l Hi £ /C*(X) the recursive sequence of compact 

subsets Hj + 2 = Fs{Hj + i,Hj), converges to A(S) with respect to the Hausdorff metric: A(<S) = lim Hj. 

j-y 00 

The natural question to make about GIFS is if they offer some new fractals. The positive answer 
is given by Mihail and Miculescu m through examples and in the recent work Strobin m for a more 
general case. We will discuss that in the end of the Section 12.11 We should mention that recently, in 
2015, Dumitru, Ioana, Sfetcu and Strobin [S] has considered many questions regarding to the extension 
of the concept of GIFS for topological contractions assuming that the family of maps is not just finite 
or countable but possibly an arbitrary family F of maps from X m to X, satisfying suitable hypothesis. 
Several results ar obtained by using code spaces (see m for details). 


1.2 GIFS with place dependent probabilities (GIFSpdp) 

In this section, we use the notation in R. Miculescu M- The set Prob(X) will always be the set of 
regular Borel probabilities on X with respect to the Borel sigma algebra induced by the metric. 

Definition 3. A generalized iteration function system with place dependent probabilities 

(GIFSpdp) is a family S of continuous functions fij : X 2 — > X, and weight functions (probabilities) 
Pj : X 2 — > [0,1] such that po(x,y ) + pi(x,y) = 1, denoted S = (X, (<fij)j= o,i> {Pj)j= o,i)- 

One special case is when the probabilities are given by a potential function u : X — > R, then 
Pj(x,y) = u{(fj{x,y)) and w((/>o(a;, y)) + u{(j)\{x, y)) = 1. We denote such case as a uniform GIF¬ 
Spdp according to Lopes and Oliveira m- 


E2 - For a GIFSpdp we assume two hypothesis on the weights: 

a) Any pfix, y) > <5 > 0 for any i = 0,1, x, y £ X; 

b) Any pfix, y) is in Lip Citdi (X 2 , [0,1]) with d + d t < 1. 

We recall that pfix,y) is Dini continuous if f 0 f^-dt < 00 for some e > 0, where Qi is the modulus 
of continuity of pi , 

I Pi{x, y) - Pi{x', y')\ < Qi ( d{{x , y), (a/, y'))) , V(x, y) £ [x !, y'). 

For instance if pi is /3-Holder ( Qfit ) = kt 13 ) or pi is /c-Lipschitz ( Qfit ) = kt) then pfix^y) is Dini 
continuous. 

Definition 4. Given S = (X, (4>j)j=o,i, (,Pj)j=o,i) we define (see R. Miculescu M)> the transference 
operator Bg : C(X, R) —> C(X 2 ,R) by 

Bs(f)(x,y)= Pj( x >y)f(((j(x,y)), 

j=o,i 

for all ( x,y ) £ X 2 . And the Markov operator Cs ■ Prob(X) x Prob(X) Prob(X) by 


[ f{f)dC s (y,v)(i) = [ B s (f)(x,y)d{g x v)(x,y), 
Jx Jx 2 

for any y,v £ Prob(X) and any continuous f : X —> R. 
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Under the hypothesis El and E2 we get, from R. Miculescu [14], Theorem 4.4, that: 

1- There is a unique ys £ Prob(X) such that Cs(ys, ys) = /is; 

2- supp(/is) = A(S), the attractor of the GIFS; 

3- For any yo,yi £ Prob(X) the sequence /ij +2 = Us (/ij, y,j+i ) converges in the Monge-Kantorovich 
distance dn □ (see [12] . Definition 2.53), to ys- 


Definition 5. TTie Hutchinson measure ys associated to a GIFSpdp is the unique solution of Cs(ys, hs) = 
hs- 

In the next sections we will consider the fixed point of another operator to get the extended Hutchin¬ 
son measure in X 2 3 . 

2 Dynamical point of view: the ergodic theorem 

A GIFS is not a typical discrete dynamical system because 4>i : X 2 —> X, is not an endomorphism 
from X 2 —> X 2 . However we can consider the dynamics of an IFS in X 2 , whose projection in the first 
coordinate, is the orbit of the GIFS. In this section, we assume the hypothesis El and E2. 

2.1 The extension of a GIFS to an IFS 

In order to analyze the orbits of a GIFS one can to embed S = (X, i n to an IFS S = 

[x 2 , (<j>i(x, y))i=o,i) where </>, : X 2 — >■ X 2 , is given by / f>i(x,y) = (y,4>i(x,y)) that is the extension 

of 5. 

Remark 6. We point out that, make a extension of a GIFS in X 2 instead X m , is not actually 
a restriction. If we consider (X,d), m,n > 2 (degree m and n maps) and a general GIFS 
as a family S of continuous functions <pj : X m —> X, denoted S — (X, (<fij)j=o...n-i), then 

his extension will be the IFS S — ^ X m , (0j(a;))j = o...n —where </>* : X m —> X m , is given 

by 4>i(x ) = (9(x), <f>i(x)), where x = (xo,x m _ 2 , x m _ 1 ) and 8 : X m —> X m ~ 1 is given 
by 8(xq, x m — 2 i x m — 1 ) = (xi,x m — 2 ). Additionally, we produce the orbits by choosing 

sequences a = (do, a-i,_) £ {0,1, ...n — 1 } *, which will make de proofs enormously hard to 

read. So, in the rest of the paper we will make the proofs for m = n = 2 making easier to 
recognize the key elements in the demonstrations 

We want to investigate the relation between the dynamics of this IFS and the properties of the GIFS 
and its Hutchinson measure. 

Definition 7. Given GIFS S = (X, (4>j)j=o,i) and a = (ao,a±,....) £ {0,1} N a fixed sequence, the orbit 
Xq,x\ € X is the sequence obtained by the nonautonomous recurrence relation xo(a) = Zo ta (xo, x\) = Xq, 
X\(a) = Zi ta {xo,x\) = X\ and Xj(a) = Zj, a (x 0 , 2 : 1 ) (or Xj = Z^ a for short) for j > 2 where -2q+i,a = 

4 , aj- 1 {Zj—l t a,Zj a '), j > 1. 

The iterations by a sequence a = (a 0 , ai,....) £ Q = (0,1} N defines the behavior of an IFS : 

{0 ao (x,y), CLQ (x,y)), •••}, 

that is c/)„, _i o • • • o <j) ai o if> ao (x, y). If we take two compact sets Hq, H 1 £ /C*(X) and a £ Hi, i = 0,1 
then 

x 2 {a) = (f ao (xo(a), Xi(a)) £ H 2 = F S (H 0 ,H 1 ), 
x 3 (a) = <f ai (xi(a), x 2 (a)) £ H 3 = F s (H 1 ,H 2 ),etc. 


J fdp - J fdv 


, for any fi, v G Prob(X). 


2 d H v) = sup 

Lips(f)< 1 

3 That is, 6(x o, a?i) = x± and 4>i(x o, xi) = ( 0(x ), cf)i(x)) = (®i, 4>i(x o, * 1 )), for i = 0,1. 
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Figure 1: Dynamics of S, the orbit of Xq,Xi by a. 


If we define Zq^(xq,Xi) = x 0 G H 0 , Z 1<a (xo,Xi) = X\ G Hi and Z m>a ( x 0 ,xi) = x m (a) £ H m , for 
m > 2, then any accumulation point of this sequence will be in H(iS), see Figure [Q Thus, this sequences 
are significant on the asymptotic behavior of the GIFS. 

Remark 8. The orbit by a GIFS S = (X, o,i) is the projection on the first coordinate of the orbit 
{(xj,Xj- |_i), j > 0} by the IFS associated S: 

(j>a 0 {x o,xi) = (xi,0 ao (xo,xi)) or x 2 = Z 2 ,a{xo,xi) = (j) ao (x 0 , Xi), 

<^ai(xi,x 2 ) = (x2,</> ai (xi,x 2 )) or x 3 = Z 3jO (x 0 ,Xi) = <f> ai {x l,x 2 ,) 

and so on, that is ^_ 1 o ■ ■ ■ o <j> ai o (j) ao (x 0 ,x i) = (xj,Xj+i), j > 1. 

In Mihail and Miculescu [IB], Ex. 4.3, we found an example of a GIFS S = (X, (</>j)j=o,i) whose 
attractor have infinite Hausdorff dimension, thus it is not an attractor of any finite Lipschitz IFS. This 

means that the GIFS theory gives us new fractal sets. In the recent work of Strobin (see Strobin ED, 

Theorem 6), the author proves that certain GIFS in X m for Icl 2 formed by generalized Matkowski 
contractions has attractors that are not attractors for any GIFS in X r with l<r<m—lin particular 
for r = 1 when the GIFS to become a classical IFS in X. However, the Proposition [T21 shows that this 
new fractals H(5) contains the projection of Al(iS), the attractor of the extension. The IFS S is not 
contractive, thus we need to consider higher powers to prove the there exists an attractor. 

Definition 9. Given an IFSTZ = {Z, (/j(z))i=o,i) we define his kth power as the IFSTZ k = (Z, (/q...i fc (z))i=o.i) 
where /q...i fc (z) = fi k o ■ ■ ■ o / ?1 (z). We said that 1Z is eventually contractive ifIZ k is contractive for some 
k > 1. Analogously, if 1Z = (Z, (fi(z))i—o,i, is a IFSpdp we define his kth-power as the 

IFSpdp lZ k where p il ...i k {z) =p ik {f ik . fh(z)) ° ■ ■ ■ °Ph(z). 

Lemma 10. The IFS S 2 is contractive, more precisely Lip(ipj) = X = m&Xj[aj + bj\ < 1 for any j = 0,1. 

In particular S has an attractor A(S) = A(S 2 ) (Or, S m is contractive and A(S) — A(S m ), if S 
has degree m.) and H n+ 1 = F§(H n ) —> A(<S) for any Hq £ /C*(X 2 ). 

Proof. We know that each <pj : X 2 — > X is in Lip aj ^ j (X 2 , X) and aj+bj < 1, that is, d(<f>j(x , y ), </>j(x', y')) < 
aj d(x,x') + bj d(y,y'). Thus, 

d((fij(x,y),<t>ij(x',y')) = d((</>i(x,y),<t>j(y, <£*(x, 2 /)))> (M x ',y'),<t>j{y >?/)))) 

= ma x[d(f>i(x,y), <t>i{x',y')), d(^j(y,^i(x,y)), <j>j(y', <f>j{x', y')))} 

< max [a* d(x,x') + bj d(y,y'), aj d(y,y') + bj d{<j>i(x,y), f>i(x',y'))\ 

< max [ai d(x, x') + bj d(y, y '), aj d(y, y') + bj {a* d{x, x') + b t d(y, y')}] 

< max [aj d{x, x') + bj d(y, y'), bjOj d{x, x') + {aj + bjbj} d{y, y')\ 

< Ad((x,y), {x',y')), 

since d((x, y), (x', y')) < maxd(x, x'), d(y , y'), where A = maxj [aj + bj] < 1. Consider the fractal operator 
Fg : )C*(X 2 ) — > K.*{X 2 ) associated to S given by F^{K) = U J= o i Ffa (-^0- 
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The operator Fg 2 is contractive because <f>i o <pj are contractions. We notice that 

Fl(K)= U F Mj (K)=F Sa (K). 

i,j= 0,1 

From Dumitru [8], Theorem 1.1, or Kunze [12] we get that there exists a unique compact set B £ /C*(X 2 ) 
such that F?(B) = Fg 2 (B) = B and H n+ \ = Fg 2 (H n ) — > B for any H 0 £ /C*(X 2 ). 

Obviously, F?(B) = B implies 0 that Fg(B) = B. The set B is the unique set with this property. We 
call A(<S) = B the attractor of S. From the fix point property we get that H n +i = Fg(H n ) —1 A(<S) for 
any Hq £ JC*(X 2 ). ■ 


In the next example we provide some computational evidence for the Proposition 1121 In these ap¬ 
proximation, we can see that A(«S) C y4(<S) 2 . 

Example 11. We consider the GIFS (f>j{x,y) = ^x + ^—p—y + for j = 0,1. A direct computation @ 
shows that <^o([0, f], [0, f]) = [0, and </>i([0, |], [0, §]) = [^,f] thus F s ([0, f] x [0, f]) = [0,^] U 
[^, |] = [0, |] = £l(5). Since the Hausdorff dimension of A(S) is one, and the Hausdorff dimension of 
proj x j4(<S) is less or equal to the dimension of A(S) that is strictly less than one it is impossible to have 
proj^S) = A(S). 

Running the correspondent chaos game with 10000 iterations we get the picture of proj a ,A(5) and 
A(S) in the Figures® and® 


23 * -*v 













Figure 2: proj a ,^4(5) C A(S) = [0, |], the attractor of S. 

Proposition 12. Let A(S) be the attractor of S, then A(<S) 2 is forward invariant with respect to the IFS 
S. Moreover, .A(<S) C ^4 (5) 2 in particular, proj x ^4(5) C ^4(5). 

Proof. Consider the fractal operator Fg : IC*(X 2 ) —1 IC*(X 2 ) associated to S. We notice that for 
K = Ki x K 2 we have F S (K i x K 2 ) = Q {K 1; K 2 ). 

i=o,i 

Since A(S) is self-similar (Fs(A(S), A(S)) = A(S)) we get 

F S (A(S) 2 )= U ^.(A(S) 2 ) = (J UA{S) 2 ) 

j=0,l 1=0,1 

= [J {(yAj(x,y))\x,y £ A(S)} 

1 = 0,1 

C |J (A(S), cfjiAiS) 2 )) 
l=o,i 

= (A(S),F s (A(S) 2 )) = A(S) 2 , 


4 Indeed, if F?(B) = B then Fg(Fj(S)) = F S (B), that is, Fj(Fg(S)) = Fg(B). Thus F$(B) = B. 
5 This computation and his consequence was pointed by one of the referees of this paper. 
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Orbit for l<n<100 




Orbit for 100<n<300 



Orbit for 300<n<1000 



Orbit for 4000<n< 10000 
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Figure 3: Building ^4(5), the attractor of S , through a random orbit. 


that is, 7l(<5>) 2 is forward invariant with respect to the IFS S. 

Taking H 0 = Hi = A(S) £ K.*(X) and Ho = H 0 x Hi = A(S) 2 £ )C*(X 2 ) the recursive sequence of 
compact subsets H n +i = Fg (H n ) C A(S) 2 thus, A(S) C A(5) 2 . ■ 


Remark 13. To study higher powers of a GIFSpdp [^| we need also consider the regularity of the weights 
in the IFSpdp S = (X 2 , (0*)*=o,i, (Pj)j= o,i) ■ In fact, if Pi{x,y) £ Lip Ci , di (X 2 , [0,1]) with Ci + d t < 1 
then, Pij(x , y) satisfy \Pij(x, y) —pij(x', y')\ < 2 d((x, y), (a/, y'). In particular each p^ is Dini continuous 
with modulus of continuity Qij{t) =2 t. 

To see this, we compute the distance 

I Pij(x,y) -Pij(x',y')\ = \Pi(x,y)pj(4>i(x,y)) - pi{x’,y')pj{<j)i{x',y'))\ 

< | p i {x,y)pj((i) i (x,y)) -p i {x,y)pj{ 4 > i (x , ,y'))\ 

+ I Pi(x,y)p j {4> i (x , ,y')) - pi{x',y')pj{(i)i(x',y'))\ 

< \pi(x,y)\ \Pj{(t>i(x,y)) ~ Pj(Mx',y'))\ 

+ \Pi(x,y) — Pi{x r , y')\ \Pj(Mx',y'))\ 

< Cjd(y,y') + djdicfiix, y),(fi(x',y')) 

+ Cid(x,x') + did(y,y') 

< Cjd(y, y') + dj[aid{x,x') + bid(y,y')} 

+ Cid(x,x') + did{y,y') 

< [ci + djai}d(x, x') + [cj + djbi + di]d(y , y') 

= q tJ d(x, x') + r i:j d(y , y), 


where qij + + djaf\ + [cj + djbi + df\ < 2. 

®We recall that the IFSpdp S 2 = ^X 2 , 

:X 2 ^X 2 

same definition holds for <S m , m > 2. 


,je{o,i} 2 


, is the 2nd-power extension of «S, where each 


(pij ■. X 2 ^ X 2 is defined by <pij{x,y) = (cpi o <pj)(x, y) = (ipj(x,y),rpi(y, <pj(x,y))), and Pij{x,y) = pi(x,y)pj(ipi(x,y)). The 












2.2 Extended GIFS and the holonomic condition 


Given a GIFSpdp S = (X, (<j>j)j= 0 , 1 , (Pj)j= 0 , 1 ) we can to embed it, in to an IFSpdp S by setting S = 
(x 2 , (Pj)j=o,i) ■ As before <j)j(x,y) = (y, y)), but we keep Pj = pj : X 2 -*• [0,1]. We 

consider the holonomic structure for an IFS introduced in Lopes and Oliveira m- Let tf : X 2 x!!-> 
X 2 x f! be the skill map 

a((x,y),w) = [4> Xo (w){x,y),(r{wfj = (y,<j) Xo ( w ){x,y),a(w)) , 

where 12 = {0,1} N , ex : S7 —>• 17 is given by a(wo,w i,...) = (wi,W 2 , ■■■) and Xfc : 12 —► {0,1} is the projection 
on the coordinate k. A cylinder is the set w o • • • w n -1 = {w £ 12 | Xo(w) = wo ,..., X n _i{w) = w n -\} C 12. 
We recall the recurrence relation xo = Zo :W (xo, x±) = x, x\ = Zi tW (xo, x±) = y and Xj = Zj tW (xo, x\) for 
j A I where Zj^.\ w — f)wj—i(Zj—\^ w ,Zj W > ), j ^ 1, so 

cr n (x,y,w) = ( Z ntW (x,y),Z n+ i tW (x,y),a n (w )) = (x n , x n+1 , cr n {w)). 

The ergodic averages evaluated on the orbits of a n allow us to define certain measures in X 2 x 12 that 
captures the behavior of the IFSpdp 


!X 2 xQ. 


h(x,y,w)drj Na(:coa;i) (x,y,w) 


1 

— h(a n (x,y,a)) 

2=0 

1 Ar " 1 

— ^ h(xi(a),x i+1 (a),cr l {a)), 
2=0 


for all h £ C(X 2 xll,l). 

Any weak limit of (lV, O!(X 0 , xl) L > 0 is a measure p a ,(xo,x{), that satisfy f X 2 xQ g ° a dy a ^ x0iXl) = 
f\ 2 xn9 ^Mo,(a;o,xi)) Vp £ C(X 2 ,M). This measures are called holonomic in Lopes and Oliveira [T51 . 
Since Prob(X 2 x 12) is compact we know that the set of those probabilities is obviously not empty. 

Definition 14. The set of holonomic probabilities with respect to the extended GIFS S = ^X 2 , (<f>i(x, y))i= 0,1 

is the set of probabilities p in X 2 x SI such that J X 2 X q 9 dp = J X 2 X ^g dp, \/g £ C(X 2 ,R). In other 
words 

[ g(0x o (w){x,y)) dp = [ g{x,y) dp, Vp£C(X 2 ,R). 

Jx 2 xn Jx 2 xn 


Lemma 15. Any holonomic measure p induces a measurable GIFSpdp S = yX 2 , (pi)i=o,ij , 

where pi(x,y) = = 0,1 are measurable weights and dp(x,y,w) = dJ^ y^w) da{x,y) is the 

disintegration of p. If p is such that pi(x,y) = J^ x ^ v ){i),i = 0,1 are continuous then S is actually a 
(continuous, finite, ...) GIFSpdp. 


Proof. By disintegration (see Dcllacherie [J) we have, for any holonomic measure p, a decomposition 
dp(x, y, w) = dJ( x>y) (w) da(x,y), 


IX 2 xCl 


h{x, y, w) dp(x, y, w) = / / h(x,y,w) dJ^^w) da{x,y), 

Jx 2 Jn 


V/i £ C(X 2 x 12, R), where the family of probability kernels Jr x ,y) is unique a almost everywhere. A 
particular case is h(x,y,w ) = g((f>Xo{w)( x :y)) for 9 € C(X' 2 ,R). Then, the holonomic condition turn in 
to 


lx 2 J a 


g{4>x 0 ( w )( x ,y)) dJ( x , y )(w) da{x, y) = 


lx 2 J a 


g(x,y) dJ( XiV )(w) da(x,y) 



g(x,y) da(x, y), 
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V <7 £ C(X 2 ,R). In other words, the marginal a(x , y) is an invariant measure for the dual of the operator 
B s(9)( x iV) = Ei=o,i J (.x, y )(j)g(fa(x,y)) ■ 


Mg £ C^X 2 ,®). 


lx 2 


B s(9)(x,y) da(x,y) = / g(x,y) da(x,y), 


lx 2 


The Markov operator associated to the extended GIFSpdp arises naturally from Lemma Upland makes 
us to consider a new class of Markov operators for a GIFSpdp. 


Definition 16. The extended operators associated to the GIFSpdp S are the usual operators to the 
IFSpdp S: 

1- Bg : C(X 2 ,R) -► C(X 2 ,R) by 


B s(9)( x ,y ) = Pj( x >y)9(4>j( x ,y)), V(x,y) g x 2 . 

i=o,i 


U- : Prob(X 2 ) —> Prob(X 2 ) by 


/ g(x,y)dCg(a)(x,y)= Bg(g)(x,y)da(x,y), 

Jx 2 Jx 2 

for any a £ Prob(X 2 ) and g : X 2 —> R. 

Definition 17. The p-operators associated to the holonomic measure dp(x,y,w) = dJ( XiV )(w) da(x,y) 
are the extended operators Bg and Cg associated to the GIFSpdp S = ^X 2 , (fa)i—o,i, (pi)i=o,i) , where 
Pi(x,y) = J(x tV ){i) = J( x ,y)({w\wo = i}). By Lemma [751 a is a fixed point, Cg(cx) = a. 

Generally speaking, the weights Pi(x, y) = J( x ,y)(i) = J( x ,y)({ w \ w o = *}) are n °t continuous. However 
it can happen if we start with the Markov process with place dependent probabilities, associated to a 
GIFSpdp (a Feller process). 

Definition 18. Given a fixed point @ a of the extended Markov operator (Cg(a) = a) associated to a 
GIFSpdp S = (X', (4>j)j= o,i, (Pj)j= o,i), one can introduce a measure in f l indexed by (x\ = y,x o = x) by 
considering his orbit x n = Z n , w (x,y) as a Markov process with place dependent probabilities. We assign 
probabilities to each cylinder wq ■ ■ ■ w n — i C ft by 


B (x,y) {wo * W n — 1 ) — Pwq (x 0 ; X\ ) • * * Pw n —i (Xn—1 ; *^n) ■ 


The product measure dp(x,y,w ) = dP^ x ^ v )(w)da(x,y), is a holonomic measure called the holonomic 
lifting (see Lopes and Oliveira m) of the measure a. Indeed, 


lx 2 x n 


g{4>x 0 (w){x,y)) dp = 


lx 2 J n 


g{f>x 0 (w) (X, y)) dP^ x y ) (w)da(x, y) 


B 


lx 2 


s ig)(x,y)da(x,y) = / g(x,y)dCgOt(x,y) 
Jx 2 


lx 2 


g{x,y) da(x,y ) = 


X 2 xQ 


g{x,y) dp. 


We recall that the IFSpdp S 2 = yX 2 , (faj). j£{01}2 , (Pij) i je{0 1}2 j , is the 2nd-power extension of S, 
where each faj : X 2 -> X 2 is defined by faj(x,y ) = (fa o fa)(x,y) = (fa(x, y), fa(y, fa(x, y))), and 
Pij(x,y) = Pi(x,y)pj(fa(x,y)). 

Lemma 19. Let S be an IFSpdp then, Cg has a unique fixed point a(ergodic with respect to the Markov 
process generated by the IFSpdp). In particular, supp(ai) = A(«S). 

7 There exists fixed points because Cg is a continuous operator and Prob(X 2 ) is compact and convex (Schauder fixed 
point theorem). 





10 


Proof. The proof uses the fact that S is eventually contractive. More precisely, we show in Lemma fTOl 
that S 2 is contractive and has an attractor A(«S). From Barnsley [I], Theorem 2.1 [f| or Kunze El, 
Theorem 2.60 and Theorem 2.63, there is a unique fixed point a £ Prob(V 2 ), for Cg 2 and supp(a) is the 
attractor of S 2 (that is equal to yf(S) by Lemma flOl) . 

Let a be any fixed point of Cg that is, CgCt = a. One can show that 

B l(g)(x,y)= Y Pi{x,y)p j (4>i(x,y))g($ i ($ j (x,y))) = Bg 2 (g)(x,y), 

i,.7=0,1 


SO, 



g{x,y) da(x, y) 


lx 2 


g{x,y)dC^a{x,y) = / g(x,y)dC\a(x,y) 

Jx 2 


lx 2 


B l(g)(x,y)da{x,y) = / Bg 2 (g)(x,y)da(x,y) 


lx 2 


lx 2 


g(x,y)dCg 2 a(x,y), 


that is, Cg 2 a = a, thus a = a. 


Definition 20. The unique fixed point of Cg in Prob{X 2 ) is the Extended Hutchinson measure for 
S. 


We can also to apply Bg in functions g(x,y) = f(Tl v (x,y)), where U y (x,y) = y is the projec¬ 
tion on the second coordinate and / £ C(X, R), obtaining Bg(g)(x,y) = E Pj{x,y)g{<i>j(x,y)) = 

3 = 0,1 


E = B s{f){x, y) that is, Bg(f oIL y )(x,y) = B s {f)(x,y), for any / £ C( X,R). 

3 = 0,1 


Lemma 21. If CgCt = a then (n x )*a: = (n y )* 0 ! that is, the marginals of the extended Hutchinson 
measure a are the same. The measure y a = (n x )*a = (n y )*a is called the projected Hutchinson 
measure. 


Proof. If we take g(x, y) = f(U x (x, y)), where n x (:r, y) = x is the projection on the first coordinate and 
/ £ C'(XjR) we have 

B s(g){x,y) = Y Pj( x iy)g(4>j( x ,y)) = 

3 = 0,1 

= E PifavKf oI[ x)(y, < f>j(x,y)) = f(y) Y Pj( x >y) = f(v) 

3 = 0,1 3 = 0,1 

that is Bg(f oH x )(x,y) = f{y), for any / £ C(X, M). Since Cg(a) = a we have 


/ f(y)da(x,y)= Bg(f oIL x )(x,y)da(x,y) = 
Jx 2 Jx 2 

= / ( f oB x){x,y)dCga(x,y)= f(x)da(x,y), 

Jx 2 Jx 2 


for any / £ Cp^K). 


8 In [T], they assume that the IFSpdp satisfy a Dini-type condition that is weaker than E2 for 2nd power extension. More 
precisely, the weights are average-contractive: YlPij hi Lip ((pi j) < 0. 

9 The subscript means the push-forward map, T*a. For T : X 2 —>■ X , the push-forward map is defined by f x fdT*a = 
fx 2 foTda. 
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2.3 The ergodic theorem for GIFS 

We address the problem of using ergodic averages of the extended GIFS to estimate the integrals 
f x 2 g{x,y)da{x,y) where a is the extended Hutchinson measure. 

Remark 22. The results will be stated for a general GIFS of degree m with n maps but the 
proofs will be made just for m = 2 avoiding the extra indexes, because there is no difference 
at all, in the reasoning. In the rest of this section, the sequence Xi(w) = Zi tW (xo, x \) for 
i = 0,1,... is obtained from ( Xq,x±) by the iteration by S. 

Our main tool is the Ergodic Lemma. 

Lemma 23. Let a be the extended Hutchinson measure. For each (xq,Xi, ...,x m -i) G X m there exists 
a measurable set H( XOtXu ... tXm _ l ) C H = (0, ...,n- 1} N such that P( Xo , Xu ..., Xm _ 1 )(£l( Xo , xu ..., Xm _ 1 )) = 1 and 
for any w G fyxo.xi.-.zm-i) 

T7 (g(x o j •••) %m— i) + ... + g(x m N —mi •••> %mN— 1 )) -> / g{x)da{x), 
tv J X 2 


V 9 eC(l m ,l). 

Proof. The proof is based on the Elton’s ergodic theorem for IFSpdp (see Elton [TTjb In order to do 
that we need to use probabilities for S 2 because it is contractive. 

Let P 2 X be the measure in f l 2 = ({0,1} x {0,1}) N indexed by (x = Xo,y = x\) 0- We built this 
measure by considering his orbit x n (w) = Z n>w (x, y) by the GIFS as a Markov process of higher order with 
place dependent probabilities (x\ = y, x 0 = x), jx 2 ,x 3 ) = f> WlWa (xp, Xi), ( x 4 , x 5 ) = f) W3W2 (x 2 ,x 3 ), ■■ ■ . 

We assign probabilities on each cylinder (wq,w\) ■ ■ ■ (w 2 n - 2 ,W 2 n-i) C H 2 (that generates the Borel 
sigma algebra of H 2 ) as follows 

P (x,y)(( W 0,Wl) ■ ■ ■ (W2n-2,w 2n -l)) = Pw 0 W! (x 0 , X \) ■ ■ • p Wn _ 1Wn (x n -l, X n ). 

The map r : H 2 —> ft given by r((wo, w±) ■ ■ ■ (w 2 n- 2 , W 2 n -i), •■•) = {wq ■ ■ ■ W 2 n-i, ■■■), is obviously a 
homeomorphism. Applying this map on cylinders with even lengths we get 


r((w 0 , Wi) • • • (w 2 n-2,W2n-l)) = Wq ■ ■ ■ W 2 n- 1 , 

and for cylinders with odd length we get 

r _ 1 (w 0 -"W2n) = (w 0 ,Wl) ■ ■ ■ (w 2 n, 0 ) U (w 0 ,W l) ■ ■ ■ [w 2 n, 1), 

so F is bi-measurable because the respective Borel sigma algebras are generated by the respective pre 
images. For each cylinder wq ■ • ■ W 2 n -1 G fl we have 

P (x,v)(( w 0 ’ w l) ' ' • (w2n-2,W 2 n -2)) = P( x>y )(w 0 ■ ■ ■ W 2n - 1 )> 

in particular T preserves measure. 

From Lemma[l9]we get that a is ergodic with respect to the S 2 . The Elton’s ergodic theorem actually 
requires an average contraction hypothesis 

Y[d(4> ij (x,y),$i j {x',y'))Pv( x ^ < r d((x,y),{x',y')), 
ij 


for some r < 1, which is a consequence of El and Lemma fTOl and that each pij (x, y) to be positive and Dini 
continuous which is a consequence of E2 and Remark [T31 Thus, we can use the Elton’s ergodic theorem 
for S 2 , that is, there exists a set f 1 2 Xq Xi ^ C f2 2 with probability one, such that for any w G f l 2 Xo s, 

(g(xo(w),x 1 (w)) + ... + g(x 2 N - 2 (w),x 2 N-i(w))) -t J ^ g(x,y)da, 

10 At this point, if the GIFS has degree m we should to consider PFJ x , the measure in = ({0,1} X ■ ■ ■ X {0,1}) N . 













12 


V g £ C(X 2 ,R). We get our result if we take I2( XOjXl ) = T(fi? a , i ^) C LI, because, 

P(x 0 ,xi){^(x 0 ,xi)) ~ ^ > (Ko>*i)(-^(^(a;o,xi))) — P(x 0 ,xi) (^(i 0 ,n)) — 


The next theorem generalizes the Elton’s Ergodic Theorem for GIFS, providing a basis for applications 
and further studies in Chaos Games Cl El El El for GIFS. This theorem generalizes the analogous result 
for IFS, Barnsley [3], Pg. 323, Theorem 4.5. 

Theorem 24. (Elton’s Ergodic Theorem for GIFS) Let a be the extended Hutchinson measure for S. 
For each (xo, x\, ..., x m -i) £ X m there exists a measurable set ^(x 0 ,xi,...,x m -i) C fi = {0,..., n — 1} N such 
that P(x 0 ,xi,...,xm—i)(^(xo,xi,...,x rn —i)') 1 and for any w £ ^l(xo,xi,...,xm—i) 

^7 (f(x K -l(w)) + ■■■ +f{xi(w)) + f{x 0 (w))) [ fdy a , 

K J X 2 

V / £ C(X, R), where y a = (II Xo )*a = (II Xl )*a = • • • = (II Xm _ 1 )*a is the projected Hutchinson measure. 


Proof. Again, we use m = 2 to simplify the writing. Consider the projections II X (x,y) = x and 
II y (x,y) = y. If we apply the Lemma[52]to g{x,y) = f(H x (x,y)) and g{x,y) = f(TL y (x,y)) we get 

((/° n y )(x 2 jv- 2 ,x 2 jv-i) + + (/ ° n x )(x 0 ,xi)) ju° n a)(a:, y)da, 

(f{x2N-i) + + fix 3 ) + f(x 1 )) — ^ J f(y)da, 

and 

-T 7 ((/on x )(x 2 jv- 2 ,a; 2 ]v-i) + ••• + (/ o n y )(x 0 , xi)) -s- / (/ o IIj / )(x, y)da, 
tv Jx 2 

(/(X2JV-2) + ... + f(x 2) + f(x 0)) —t j f(x)da. 

Adding this two limits we get 

2 ^ (f{x2N- 1 ) + f[x 2 N-2) + ... + fix 1 ) + fix 0 )) J /(x) + /(y)do, 

thus 

(/(aw) + + /(®i) + fi x o)) ^ J fix ) + /(y)da(x, y), 

if K = 2IV. On the other hand, if AT = 27V + 1 we use the fact that goes to zero uniformly, what 
guarantees the result for any K. Since (II x )*a = y a = (II y )*a, the limit above will be 

if{x K -l{w)) + ... + /(xi(u>)) + f(x 0 (w))) -> J f{x)dy a {x). 


Remark 25. A natural question is if a = yx y. The answer is not in general, but it can hapens in some 
cases. One can easily test this possibility using functions gix,y) = a(x)6(y). Indeed, even for a standard 
GIFSpdp in X = [0,1], <j) 0 {x,y) = \x + \y and fa (x,y) = \x + \y + \ and Pjix,y) = |, j = 0,1 we 
for a(x) = x and a generic b, an impossible equation appears, if we suppose a = y x y. In the 
Example\W\ of the Section [70 we found a = ys x ys- 


3 Applications and Examples 

In this section we use the ergodic theorem to get several consequences and applications to related fields. 
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3.1 Chaos Game: a random iteration algorithm for GIFS 

In the 80’s, M. F. Barnsley [HU has introduced the idea of Chaos Game or Random Iteration Algorithm 
as a tool for drawing fractals appearing as attractors of IFS. However, in the last few years, the term 
Chaos Game means every iteration (sets, multifunction, point) using a random choice of maps. The 
algorithm is described as follows (see Barnsley [3], Chapter III, or Kunze [12] , Chapters 2 and 6 for more 
details): 

Let 1Z be an attractive IFS on X, that is a family of functions <pj : X —> X, and we introduce 
weight functions ( probabilities) pj : X —>• [0,1] such that po{x) + p\(x) = 1 producing an IFSpdp, denoted 
IZ = (X, o,i, (Pj)j= o,i)■ If A is the attractor of S then, for each xo £ X 

A= lim {x n (a)} n >k, 

k—> oo 

for almost every a £ fl = {0,1} N with respect to P Xo Where the orbit x n (a ) is obtained from Xq taking 
Xn+i = 4>a n { x n) and a n is chosen with probability Prob(a ra = j) = pj(x n ). 

We will show that the chaos game, in the sense of approximate the picture of the attractor by the 
closure of a random orbit, works for a extended GIFS. However, as we can see from Example [TT] that 
satisfies the hypothesis of Theorem [26] we are not able to draw the attractor of GIFS by projecting the 
attractor of his extension. 

Theorem 26. (Chaos game for extended GIFS) Let S be a GIFS satisfying El and E2 hypothesis. If 
A(5) is the attractor of his extension then, for any fixed xq,x\, ...x m -i € X we have 

^(t^) — lim { (x n (^), • ■ •, x n+m— 1 (a) )}n>A: ? 
k—>oo 

for almost alia £ Q with respect to the probability P Xo , x where {x n (a)} n is the orbit of(x o, x\, ...x m _i). 
In particular pro) X A(S) = lim { x n (a)} n >k ■ 

k—> oo 


Proof. We assume m = 2. To prove that A(<S) = lim^oo {(x„(a),x n+ i(a))}n>fc, we take 2 : £ A(<S) = 
suppa then there is a neighborhood V of z such that a{V) > 0. By Lemma H3H there exists a measurable 
set C H such that P( X0 ,x 1 )(^{x 0 ,x 1 )) = land for any a € ^l( Xo , Xl ), \ (xj(a), x j+ i(a)) £f}-> 

a(V). Reducing the size of V we get that z £ lim^oo {(xj(a),Xj+i(a))} n >k thus A(S) C lirn^oo {(xj(a), Xj+i(a))} n >k- 
On the other hand, A(5) is forward invariant. If we have (xo, Xi) £ A(£>) then A(<S) A limfe^oo {(xj(a), Xj+i(a))}„>fc, 
what proves the equality because the closure of the orbit does not depends on the first point, only on 
the sequence a £ H. The second part proj a ,A( l S) = lim^oo {x n (a)} n >k follows in the same fashion using 
Theorem 1241 or we can use the first part of this proof and the fact the projection is continuous. ■ 


3.2 Nonautonomous Dynamical Systems 


A well known example of Discrete Nonautonomous Dynamical System (see Potzsche HU) is the dynamics 
generated by a finite difference equation(FDE), defined by a nonautonomous recurrence relation of order 
m > 2, nominally 

Xj+m — f(.Xj-\-Tn— 1 5 Xj^-jyi —2 j • *■ j Xj , CLj ), J A 0, 

where ao,ai,... £ I. The control sequence (fflj)j=o,i,... represents some seasonal interference acting on 
each iteration by changing the standard recursiveness on a FDE. It can be modeled by a GIFS S = 
C X m , o...n-i) where <f)j : X m X is given by 


(VOi 2/1 ) Urn— 1) — /(Z/0i J/li 2/m—1; Ojj). 


The control set I could be finite or not. In this case, the orbit of the GIFS from the point (co,..., c m _ 2 , c m _i) 
is equal to the orbit of the FDE with initial conditions: 


Xo 


— f 1 •••? %j 1 Cl (z I 

= Cq 


•Em— 1 1 — Cjji— 1 


11 The limit is taken with respect to the Hausdorff metric. 
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If the associated GIFS satisfy the hypothesis El and E2, we can apply our theory to study the asymptotic 
behavior and the limit sets of these nonautonomous dynamical systems. 

Example 27. A bad example is X = R, noncompact, and the not contractive GIFS given by 4>-\(x,y) = 
x — y and <j>i(x,y) = x + y or 4>j{x,y) = x + jy associated to the FDE 

{ Xj -\-2 — ./; Xj , of) — .r'j'+i F CLjXj 

Xq , Co 

Xu = C 1 


Here,I = {—1,1} and 11 = {—1,1} N . In this case the theory does not work because there is no global 
attractor. However A= {(0,0)} is an invariant set. 

Example 28. A good example is the second order FDE in X = [0,1] given by 


x 0 +2 = \x j+1 + \xj F aj =j£ {0,1} 

x 0 , = c 0 G [0,1] 

xi, = ci G [0,1] 


Associated to this FDE we have S = (X,(<f>j)j = o,i) an GIFS where (f>j{x,y) = + Both functions 

<pj are Lipi i ([0, l] 2 , [0,1]) (this example appears in R. Miculescu w- 

The typical limit set of the orbits ( x n , x n +i), of this FDE is the attractor of the corresponding extended 
GIFS. In this case we can show explicitly a formula for those points: 


Xn+2 


On+1 

4n+l 




4_n +1 


Xi 


4 n 2 


a o 


4° 2 


n > 0, 


*715 ' v _ 


where x 0 = Z 0 , a {xo,xi), x\ = Z l!a (xo,Xi) and Xj = Zp a (x 0 ,Xi) satisfies x j+1 = 4> aj -i{ x j-i,Xj), j > 1. 
And, the numbers bg, b\,... are given by the generating function G(w) = x —a = 1 + u> + 5 w 2 F 9 w 3 + 
29 w 4 + 65 w 5 + 181 w 6 F 441 w 7 + 1165 w 8 F O ( w 9 ) . Solving 1 — w — 4mj 2 = 0 with respect to z we get 
b n = -^= [A^" n_1 — A)~ ra_1 ] , n > 2. Finally, the coefficients of xq and x\ are of the type 


bn— 1 

22n+2 


1 

2 2n + 2 Vl7 


A - l ™ -1 ) -1 _ y - l 71 - 1 )- 1 


Ao 


Ar 1 


4y/rf L( 4 A 2 ) n (4A 1 ) n _ 


because ||4A,;|| > 1, i = l,2(4Ai « —2.56... and 4 A 2 ~ 1.56...). 

By the Chaos Game Theorem, the closure in X of the orbit defined by a = (ao,ai,...) is proj a ,j4( l S): 


proj^S) 



= lim cl 

k—> 00 



b j 

41 a ™~ 2 ~i 


for a random a G 12. The picture of A(S) is given by the Figure ^ and is drawn using the Chaos game 
theorem. 

Employing the ergodic theorem we can also estimate the frequency of visitations of the solutions of 
the FDE on a subset of the phase space. We choose an initial distribution po(x,y) and pi(x,y) satisfying 
hypothesis E2 (such as, po(x,y) = 4 and pi(x,y) = |{. From Theorem m if p a is the Hutchinson 
measure for S and B C [0,1] there exists a measurable set 12( COjCl ) C Q such that P( C0;Cl )(12( C0>Cl )) = 1 
and for any a G ft( Co>Cl ) 

{o<j<iv —1 | Xj (&) G y 

That is the typical average of visits of the orbits on the set B, with respect to probabilities that we choose. 


3.3 An example from Thermodynamic Formalism 

In Thermodynamic Formalism (see Bowen [5] or Walters [22] for more details), the Ruelle theorem gives 
Gibbs measures for a potential with respect to a map. 
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Figure 4: suppa = H(<S), the attractor of S. The picture of the attractor of extension is obtained by 
4000 iterations of a random orbit. 


Theorem 29. (Ruelle) Given, (X,d) a compact and connected metric space, T : X —> X an expanding 
map n to 1, A : X —> R a Holder potential, there exists a positive number X, a function h and a measure 
v such that PAh = A h and P\v = \v , where, Pa{/){x) = Ylry-x e" 4 *- 1 '■’/(?/) ® s ^e Ruelle operator and his 
dual Pj^, act in probabilities. Moreover, p = hi/ called the Gibbs measure, is T-invariant. In particular 
the Gibbs property implies that p positive in open sets. 

Usually, such dynamical systems are identified as a uniform IFSpdp 
U = {X, (pi(x) = 

where the inverse branches of T are the contractive maps tj, that is, T o Ti(x) = id(x ), i = 1, 2, n. In 
Jorgensen [10] this identification is called the endomorphism case. If the potential is normalized, that is, 
A = 1 then h = 1 and the Gibbs measure p = pn is the Hutchinson measure of the IFSpdp because 

p A(f)(x) = eMv) f(y) = Y Pi( x )f(n(x)), 

Ty—x i—l..n 


is the fractal operator associated to 1Z. 

The easier case where this happens is for X = [0,1], T(x ) = 2x mod 1 and A(x) = In T p x ^ = — In 2. 

In this case the inverse branches are Tj(x) = \x + i = 0,1 and pt{x) = e -ln2 = It is well know, 
that the Gibbs measure is the Lebesgue measure dx in the interval (see Conze [6] for a detailed study of 
the operator Pa in this case). 

The next example shows that ^4(5) = Al(<S) 2 can happen. We consider a GIFS that “duplicates” the 
behavior of the IFS 1Z built from T(x) = 2x mod 1. 

Example 30. The GIFSpdp S in [0, 1] given by </>i(x,y) = hx + i = 0, 1 and Pi(x,y ) = 5 , i = 0,1, 
satisfy: 

a) ds = dx and A(S) = [ 0 , 1 ]; 

b) a = dxdy and T(5) = [0, l] 2 . 
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Figure 5: Chaos game, for 40000 iterations, showing that A(S) = A(<S) 2 . 


Indeed, 


f(x)Cs{dx , dx) = 


J j Bs(f)(x, y)dx dy 

J J Po{x,y)f(<l>o(x,y))+Pi(x,y)f{<j>i{x,y))dxdy 

J S\ ! {r + t} + \ ! {\ x+ ^) ixiv 

j J P-\ n 2 {f)(x)dxdy = J dy J f(x)Pl ln2 (dx) 

[ f(x)dx, 


thus Ls(dx,dx) = dx, what means that the Hutchinson measure ys is equal to the Lebesgue measure dx. 
Since A(<S) = supp ys we get A(S) = [0,1]. 

In order to prove the second claim, we consider the extended operator 


-II 

-II 

-II 

-I! 

-II 


B s(9)i.x,y)dxdy 


Po(,x, y)g ($o (x,y)) +p 1 {x,y)g (^i(x,y)) dxdy 


1 

f 1 

0\ 

1 

( 1 

i\ 

2 51 

[ y '2 X + 

2 

1 + 2 9{ 

\y^ x + 

2) 


P-\n 2 (g(y, -)){x)dxdy = J J g(y,x)Pl ln2 (dx)dy 
g(y,x)dxdy = J J g(x,y)dxdy, 


thus Cg(dxdy) = dxdy, what means that the extended Hutchinson measure a is equal to the Lebesgue 
measure dxdy in [0, l] 2 . Since A(S) = suppa we get ^4(<S) = [0, l] 2 . 
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